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Abstract 

We uncover a surprising correspondence between a non-perturbative formulation of 
three-dimensional Lorentzian quantum gravity and a hermitian two-matrix model 
with ABAB-interaction. The gravitational transfer matrix can be expressed as the 
logarithm of a two-matrix integral, and we deduce from the known structure of the 
latter that the model has two phases. In the phase of weak gravity, well-defined 
two-dimensional universes propagate in proper time, whereas in the strong-coupling 
phase the spatial hypersurfaces disintegrate into many components connected by 
wormholes. 
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1 Introduction 



A major task of modem theoretical physics is to unite quantum mechanics with the 
theory of gravity and to understand "quantum geometry" . In four dimensions, this is 
proving a difficult task and there is no general consensus on which direction to take. 
The so-called "string community", originating from quantum field theory, claims 
that M-theory provides the only viable road to unifying quantum mechanics and 
gravity, while people coming from the theory of relativity tend to favour approaches 
based on canonical quantization. So far neither of these approaches has given us a 
detailed understanding of the microscopic quantum geometry of the real space-time. 

The situation is considerably better in dimension d < 4. Although lower- 
dimensional models do not possess propagating gravitational degrees of freedom, 
their geometries are still subject to quantum fiuctuations, and the quantum the- 
ories are non-trivial. One may use conventional quantum field-theoretic methods 
to investigate the coupling between matter and gravity, and to define and calcu- 
late diffeomorphism-invariant correlation functions and the dynamically generated 
fractal dimension of quantum space-time. Of course, we have no way of knowing 
how relevant the study of the lower-dimensional theories will be for the eventual 
theory of four-dimensional quantum gravity, but it is clear that there are structural 
similarities and certainly some of the same questions can be asked. 

A detailed and explicit analysis exists in two space-time dimensions. One can 
calculate the anomalous dimensions acquired by matter fields when coupled to 2d 
Euclidean gravity 0. It is also understood that a typical two-dimensional Eu- 
clidean geometry contributing to the gravitational path integral has a fractal struc- 
ture with Hausdorff dimension four, much in the same way as a typical path in the 
path integral for a particle is also fractal, with Hausdorff dimension two. 

The fractal structure is best understood by introducing a geodesic "time" on the 
two-geometries. Euclidean 2d gravity is characterized by the fact that an infinite 
number of baby universes branches off the one- dimensional spatial slice as it evolves 
in this "time" [@, i|, H]- While such a process is unavoidable in two-dimensional 
gravity models coming from string theories (describing non-critical strings), there 
is nothing in a theory of quantum gravity which demands that space should be 
allowed to split into disconnected parts. In the context of canonical gravity this 
would require a "third quantization" to enable the destruction and creation of baby 
universes, a possibility that is not usually considered. 

Not allowing for the creation of baby universes leads to a new two-dimensional 
quantum theory, called Lorentzian quantum gravity. The name derives from the fact 
that the sum over geometries in the path integral includes only a particular subclass 
of Euclidean geometries which are obtained through a Wick rotation from a set of 
Lorentzian space-times with a well-defined causal and globally hyperbolic structure 
. Interestingly, the resulting quantum theory is different from the Euclidean one. 
For instance, the fractal dimension of a typical geometry is two and not four, and 
the coupling of matter and geometry creates no anomalous scaling dimensions for 
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the matter fields [0, The detailed relation between the geometries of the two 
models is well understood 

Very important in the study of two-dimensional quantum gravity have been 
methods and concepts from statistical physics JTO], |]. The regularization in both 
the Euclidean and Lorentzian case uses so-called dynamical triangulations. In this 
approach, geometries are created by gluing together large numbers of identical trian- 
gular building blocks. The geodesic edge length of the triangles is a measure of the 
fineness of the simplicial lattice and defines a diffeomorphism-invariant cutoff of the 
theory. The action of such a piecewise linear two-geometry is calculated by Regge's 
prescription, and the state sum over geometries can in many cases be performed 
explicitly. The scaling limit (i.e. taking the number of triangles to infinity), defines 
the continuum limit of these models. 

Since in both Euclidean and Lorentzian 2d dynamical triangulations we have a 
concept of time ("geodesic" time and "proper" time respectively), a transfer matrix 
can be introduced which describes the evolution (the transition amplitude) between 
the spatial configurations at time t and t + a, where a is a discrete lattice spacing. 
The notion of a transfer matrix is familiar from quantum field theories on fixed 
lattices. It allows us to extract the continuum Hamiltonian of the system in the 
limit as a — > according to 

(0(a:)|T|0'(x')) = (0(a:)|e-<^^|0'(x')) - (0(a:)|(i - ai^ + O(a2))|0'(x')). (1) 

This strategy has been applied successfully in Euclidean and Lorentzian two-dimen- 
sional quantum gravity. 

While all of this works beautifully in dimension two, where it has given us a num- 
ber of powerful analytical tools, the situation is quite different when one tries to use 
the method of dynamical triangulations to obtain a theory of Euclidean quantum 
gravity in higher dimensions. Such models have been investigated mainly through 
numerical simulations |TT| , |12|, and the results have so far been disappointing: no 



interesting continuum limits seem to exist [|T^.0 In fact, this was one of the main 
motivations for constructing alternative Lorentzian models of dynamical triangula- 
tions, which we have already mentioned in the two-dimensional context. They have 
recently been shown to exist as well-defined regularized models of quantum gravity 



also in three and four space-time dimensions |jT5|, (see also [|T^ for a review of 
discrete Lorentzian gravity). 

An investigation of the continuum properties of the three-dimensional Lorentzian 
gravity model has already begun. Computer simulations show that it avoids some 
of the problems of Euclidean simplicial quantum gravity and most likely has a con- 



tinuum limit ||18[, thus fulfilling some of the hopes raised by the two-dimensional 
Lorentzian model. As is well-known from previous attempts, analytic tools are hard 
to come by in statistical models of quantum geometry m d > 2. In particular. 



^This seems to be part of a general pattern, since also other non-perturbative discrete approaches 
to four-dimensional Euclidean quantum gravity have had little success, see for a review. 
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the matrix-model methods that proved so powerful in two dimensions have not yet 
been made into a useful calculational tool in higher-dimensional quantum gravity. 
The purpose of this article is to demonstrate that in three-dimensional Lorentzian 
quantum gravity, such analytic matrix-model techniques can indeed be employed. 

The remainder of this article is organized as follows. In Sec. ^ we describe 
the quantum gravity model in terms of simplicial geometries and define its trans- 
fer matrix. At the discretized level, it is a simple variant of the three-dimensional 



Lorentzian model introduced in [|l^, |16| and studied by Monte Carlo simulations 
in [p!8| , p!9|] . Motivated by some well-known properties of (2-|-l)-dimensional quan- 



tum gravity, we perform an integration over all but one of the spatial geometric 
degrees of freedom in Sec. ^ In the following Sec. § we remind the reader of how 
to obtain the continuum Hamiltonian from the resulting transfer matrix. The cor- 
respondence of 3d Lorentzian quantum gravity with the already partially solved 
two-matrix model pO[ with y4i?y4i?-interaction is established in Sec. |^. We rein- 



terpret the phase structure of this matrix model in terms of geometry in Sec. g. 
This also involves a discussion of other, closely related matrix models with so-called 
touching-interactions. We briefly comment on the status of the full and as yet 
unsolved ABAB-matiix model in Sec. |^. In Sec. ^ we describe how taking the con- 
tinuum limit in the ABAB-model fits in with our previous considerations of this 
limit in three-dimensional gravity. We end with an outlook in Sec. ^. The appendix 
contains a derivation of the gravitational action in terms of the 3d building blocks 
used in this article. 



2 Lorentzian 3d gravity from pyramids and tetra- 
hedra 

The motivation for constructing non-perturbative gravitational path integrals of 
Lorentzian geometries and the general properties of the dynamically triangulated 



model in three dimensions were described in ||15| , p!q] . For our present purposes it 
is convenient to consider a slightly modified regularization which can be related to 
a quartic matrix model. (The cubic matrix model that would correspond to our 
original model which uses only tetrahedra has not yet been solved.) 

In the regularized model, (proper) time t is discretized into integer lattice steps of 
unit one. The spatial slices at t = 0, 1, . . . are piecewise linear manifolds of spherical 
topology, constructed by gluing together fiat squares with edge length Ig = a (rather 



than the equilateral triangles with edge length Is = a of references JTS] , 



The geometry of a spatial slice is uniquely fixed by this length assignment and by 
its connectivity matrix, specifying which pairs of squares are glued together along 
a common edge. Two spatial slices at t and t + 1 form the space-like boundaries of 
a three-dimensional piecewise-linear manifold "sandwich" that lies in between the 
slices and has topology S"^ x [0,1]. The fundamental three-dimensional building 
blocks used for "filling in" are regular pyramids with square base and tetrahedra 
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Figure 1: Pyramids and tetrahedra can be used to discretize 3d Lorentzian space- 
times. We show the three types of fundamental building blocks and their location 
with respect to the spatial hypersurfaces of constant integer-t. 

(see Fig. |1|). The pyramids have either their base in the t-plane and their tip in 
the t + l-plane - in which case we call them (4,l)-pyramids - or vice versa for the 
(l,4)-pyramids. In addition, each pyramid has four time-like edges of equal squared 
length connecting the neighbouring slice^. Two pyramids of equal orientation 
can be glued together along a time-like triangle. Since (4,l)-pyramids with base at 
t and (l,4)-pyramids with base at t + 1 cannot share a triangular face, we need an 
additional type of building block, namely, a (2,2)-tetrahedron with one space-like 
edge (of length Ig) each in the t- and t-|-l-planes (Fig. Its remaining four time-like 
edges have again squared length l^, so that its faces can be glued to the pyramids 
of both types. 

We will in the following use the term "quadrangulation" to denote a piecewise- 
flat geometry made of (4,1)-, (1,4)- and (2,2)-building blocks, as well as their two- 
dimensional spatial sections. (Note that a 3d quadrangulation may be thought of 
as a particular kind of a three-dimensional trzangulation, obtained by cutting all 
pyramids into pairs of tetrahedra.) Starting from a ^^-quadrangulation at time 
t = 0, we can by successive gluing build up a three-dimensional space-time of length 
t in the time direction, consisting of t + 1 spatial spheres, and t "sandwiches" in 
between. 

Although our choice of allowed discretized space-times is of course motivated 
by the causal structure associated with the physical Lorentzian signature, we will 
from now on do all calculations for the already Wick-rotated Euclidean geometries. 
Without loss of generality, we set = = /5a^, for > 1/2 (because of triangle 
inequalities). Next, we must compute the Boltzmann weight e~^ associated with a 
given three-geometry. The Euclidean Einstein action in the continuum is given by 

5=-* = -tttV / d^^V^(^(a:) - 2A) - / d'x^/hix)K{x), (2) 

■^Our language here is "Lorentzian" in line with our general philosophy [ p^ p^ , where a Wick 
rotation corresponds to a sign flip of l^. 
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where h{x) is the induced metric and K{x) the trace of the extrinsic curvature on 
the boundary dM of the manifold M, and and A are the gravitational and 
cosmological coupling constants. As usual for 3d piecewise linear manifolds, the 
curvature is concentrated on the one-dimensional edges or links, and proportional 
to the deficit angle under rotation around a link in a plane perpendicular to it. 
Using the standard Regge prescription for computing the total scalar curvature of a 
simplicial manifold, the discrete counterpart of is derived in the appendix, where 
we also discuss the inverse Wick rotation of the action. 

Let us adopt the following notation: spatial quadrangulations at integer-t are 
called T{t) and a sandwich geometry Q{T{t),T{t + l)). One piece of information 
contained in the data characterizing the geometry Q{T{t),T{t + l)) in the interval 

+ is the set of three numbers A^4i(t), A^i4(t) and A^22(^) of building blocks of 
the three types. For a given geometry Q, the action depends only on these three 
bulk variables, 

S[giTit),T{t+l))] = co-fc(Ar4i(t)+Ari4(t)-Ar22(t)) +A(Ar4i(t)+Ari4(t)+iAr22(t)). 

(3) 

In (|^), A^4i(t) + A^i4(t) — A^22(^) is proportional to the integrated scalar curvature 
(including the extrinsic curvature terms for the boundaries) between t and t+1, while 
N^iit) + Niiit) + ^N22(t) is proportional to the three- volume of this piece of space- 
time. The dimensionless coupling constants k and A are proportional to the bare 
inverse gravitational coupling constant I/Gn and the bare cosmological coupling 
constant A. The explicit form of the coupling constants as functions of (3, Gn and 
A can be found in the appendix. 

When we start stacking up sandwich geometries, the action @ is by construction 
additive, so that the total action for a space-time extending t steps in the time- 
direction becomes 

S[k, A] = Cot - ~k{N,, + Nu - N22) + A(iV4i + Nu + ^iV22), (4) 

where A^4i, A''i4 and A^22 denote now the total numbers of (4,1)-, (1,4)- and (2,2)- 
building blocks. 

We can now define the transfer matrix associated with a unit step of proper time 

asQ 

{T{t+l)\f\T{t)) = Yl ^e-^[^(^«'^(«)W, (5) 
g{r{t),r(t+i)) ^ 

^ Following [p^ , we introduce quantum states |T) at fixed i, labelled by inequivalent spatial 
geometries and normalized according to 

(Ti|r2) = -^5r,.r„ ^Cr|r)(T| = i, 
where Cr is the symmetry factor of the quadrangulation T. 
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where the summation is over all distinct three-geometries Q{T{t),T{t + l)) whose 
boundary geometries are T{t) and T{t+1). Cg denotes the symmetry factor of the 
quadrangulation Q, i.e. the order of the automorphism group of Q. 

Let Ti = T{t) contain A^4i(T(t)) squares and T2 = T(t+1) contain Ni4^{T{t+l)) 
squares. Substituting the sandwich action into (j|), we can write 

22 

(TalflTi) =e-^°e-M(^^^(^^)+^^*(^^» J2 Ar(Ti, Ts, iVsa) e-(5^+^)^22^ (g) 

where the sum is over all values N22 which can occur in sandwich geometries with 
boundary Ti U T2. The combinatorial factor M(7i,T2, N22) counts the number of 
distinct three-geometries (including the symmetry factor weights) in for a 

fixed number of (2,2)-tetrahedra. 

The Euclideanized amplitude for propagating a spatial geometry Ti = T(0) at 
proper time to a later geometry T2 = T{t) at proper time t is obtained by a t-fold 
iteration of the transfer matrix, 

G(ri,r2,t) = (T2|f*|ri), (7) 

and satisfies the completeness relation 

(T2|f*|Ti) = ^(T2|f*i|T)Cr(T|r*^|Ti), (8) 
r 

for any split t = ti+t2 of the total time interval. 



3 Integrating out geometries 

Let us first recall the situation in classical 3d gravity, on space-times with compact 
spatial slices S^^^ of genus g. The degrees of freedom of the theory associated with 
any spatial slice are the geometries (i.e. the spatial metrics gij modulo spatial dif- 
feomorphisms) , the elements of superspace. A priori, they are genuine field degrees 
of freedom: for example, each metric can be decomposed uniquely (up to a diffeo- 
morphism) into a constant-curvature metric cjij and a conformal factor according to 
gij{x) = e'^^^^^gij{x). However, a canonical analysis reveals that the conformal factor 
\{x) is not a dynamical field degree of freedom, but is completely determined by 
solving the constraints. What remains is a finite number of so-called Teichmiiller 
parameters (none for = 0, 2 for g = 1 and 6(7 — 6 for g > 1), coordinatizing the 
space of constant-curvature metrics for a given compact spatial manifold of genus 
g. This implies that for the spherical case with E*^°) = S'^ not even a finite number 
of (classical) degrees of freedom is left after the conformal factor has been fixedQ. 

^We mean here "bulk" degrees of freedom; depending on the choice of boundary conditions on 
the initial and final spatial slices, surface degrees of freedom may be "liberated". 
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This raises the question of how these classical properties are reflected in a gra- 
vitational path-integral approach of the type we are considering. At the discretized 
level, we can in principle compute the quantum amplitude between two arbitrary 
spatial geometries (either at a fixed proper-time distance t, or for an arbitrary dis- 
tance, obtained by summing the discrete propagator (|^ over all positive integers t). 
In general, this amplitude will not vanish, since there are always many 3d quadran- 
gulations interpolating between two given boundaries Ti and T2. It implies that a 
generic path (i.e. a three-geometry) contributing to the path integral does not obey 
the classical constraints, which is not particularly surprising. However, what one is 
really interested in is the behaviour of these amplitudes in the continuum limit. 

Previous investigations of the non-perturbative gravitational path integral^ are 
suggestive of what may be happening in three dimensions. The continuum analysis of 
[p3| found that - subject to a number of plausible conditions - the kinetic term of the 



conformal factor A is cancelled in the non-perturbative path integral by a Faddeev- 
Popov term in the measure. This means that there is no conformal kinetic term 
in the effective action, and therefore that X{x) is not a propagating field degree of 
freedom. This scenario is corroborated by our numerical simulations of 3d Lorentzian 
gravity []TB|, |TP[ which did not show any evidence of the conformal divergence closely 



associated with the kinetic term for the conformal factor. 

Motivated by these considerations, we will introduce a vastly reduced set of 
(Hilbert space) states, by summing over all geometric degrees of freedom of a spa- 
tial slice of a given two-volume A^. Writing iT/v) for a state corresponding to a 
quadrangulation with squares, we define for a fixed spatial topology the state 

\N):=^^=J2CrM- (9) 



The normalization factor in front of the sum ensures the orthonormality of the 
"area states", (A^|A^') = 5n,n'- The discrete area A^ may be thought of as a global 
conformal degree of freedom. We do not integrate over A^ in order to keep control 
over the continuum limit and to be able to compare our results with those of similar 
continuum approaches (see, for example, ^)- Using the new states (|^), we now 



make the conjecture that 

(Ari4|f|r;v4i)-(iVi4|T|r4^)^0 for Ari4,Ar4i^oo. (10) 

Expressed in words, this means that for large areas A^4i and A^i4 the expectation value 
(A'i4|T|7/v4i) does not depend on which "representative" Tfq^^ is chosen from the set 
of 2d quadrangulations with A'41 squares. The conjecture should be understood in 
a probabilistic sense: it says that the number of states where (|10D is not satisfied 
should have a slower growth as a function of A^ than the total number of states |T/v) 

^What is relevant for our current purposes is always a configuration space path integral, and 
not one in phase space. 
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(which grows exponentially with for a given spatial topology). Heuristically one 
can view ([TUp as expressing that in the large- limit the matrix elements {T]^\T\Tl^i) 
depend not on the two quadrangulations separately, but only on a suitably defined 
"distance" between T/v and 7^,, similar to the way the integral / dxf{x—y) over the 
real line is independent of y. In principle it is a combinatorial problem to show that 
(p!0D is valid, but we have not yet produced such a proof. 

An important consequence of property ([T0|) is that the completeness relation (j^) 
continues to hold for the area states |A^), 

(iV2|f*^+*2|A^i) = ^(A/'2|f*2|A^)(A^|f*i|A^i), (11) 

N 

where we have again assumed that A^i, A'2 and A^ are all large. We therefore still 
have a transfer matrix formalism, but with the transfer matrix T acting only on the 
subspace spanned by the linear combinations {|A^)} of the original Hilbert space. 
Obviously, the task of diagonalizing the transfer matrix on this reduced space is 
considerably simplified compared with the original problem. 

From now on, we will focus our interest on solving the combinatorics of a single 
sandwich geometry. We can rewrite relation (|^) in an obvious notation as 

/\rmax 

(Ari4|f |Ar4i) = e-^«^(^-^)(^«+^") ^ AfiN^i, N^, N22) e-'^'^'^+''^^^\ (12) 
where 

Af{N^uNM,N22)= J2 ^{'^N,,,Ti,,,,N22) (13) 

denotes the total number of quadrangulations of the space-time between t and t+l, 
including sums over the connectivities of the spatial boundary geometries at t and 
t + 1. One may think of (0) as describing the combinatorics of quadrangulating a 
sandwich geometry with free boundary conditions, except for the areas of the two 
boundaries which are kept fixed. 

In order to simplify the combinatorics further, we introduce boundary cosmolog- 
ical constants Aj and Af associated with the initial and final boundaries at t and 
t+1. They do not have an immediate physical interpretation and should simply be 
thought of as convenient book-keeping devices that will be set to zero at the end|]. 
They give rise to an additional term 

AS = KalN^i + A/a^A^i4 = ZiN^i + -Z/A^u (14) 

in the action (^, which will allow us to introduce an asymmetry in the coupling 
constants multiplying the two areas A'41 and A'^14. We can use the dimensionless 

^This is analogous to the introduction of external sources into a quantum field-theoretic path 
integral in order to obtain Green's functions through functional differentiation. 
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boundary cosmological terms Zi and zj to obtain the Laplace transform of the trans- 
fer matrix with respect to both the initial and final areas, Ni = A^4i and Nf = N14, 
namely, 

{zf\f\z,) = J2 e-^^^-^^^^(iV;|f |iV,). (15) 



f 



The matrix elements (A^j|T|A^j) are calculated for Zi = Zf = 0. From a combinatorial 
point of view, (A^/|T|A^j) counts the three-geometries Q for given numbers Ni and 
Nf, each with relative weight e"'^^'^"'"'^)^^^^ g^^d {zf\T\zi) plays the role of a generating 
function for these numbers. It is usually much easier in combinatorial problems to 
calculate the generating function rather than the actual numbers. This is well illus- 
trated by both Lorentzian and Euclidean simplicial quantum gravity in dimension 
two, where boundary cosmological constants are introduced in an analogous man- 
ner, and where the associated generating functional greatly simplifies the counting 
of geometries. 

The reason for the simplification is the fact that by going to {zf\T\zi) we have 
achieved totally free boundary conditions, since the constraints of fixed areas iVj 
and Nf have been lifted. On the other hand, no information has been lost, since 
we can in principle always rederive {Nf\T\Ni) by an inverse Laplace transformation 
from {zf\T\zi). Given the one-step propagator (A^j|T|A^j), we can finally obtain the 
propagator G{Ni, Nf,t) for arbitrary times by iterating according to ([111) . 

4 Extracting the Hamiltonian 

Once the matrix elements {Nf\T\Ni) are known, one may try to extract the con- 
tinuum Hamiltonian operator H of the system by expanding them in the lattice 
spacing a and then taking a ^ 0, 

{Nf\f\N,) = {Nf\e-''^\Ni) = {Nf\(i - aH + 0(0^)) \N,). (16) 

In this way one obtains the quantum Hamiltonian in the 'W-representation". It is 
also possible (and usually easier) to extract H from the Laplace transform of the one- 
step propagator, {zf\T\zi) (yielding the Hamiltonian in the "dual" z- representation). 

Let us illustrate this by a concrete calculation in 2d Lorentzian gravity, where the 
space-time has topology x [0, 1]. In this case, the generating functional {zf\T\zi) 
is known explicitly^ 0, 

{zf\T\z,) = log ^ . ■ (17) 



^In line with our construction in dimension three, we consider here the symmetric propagator 
of with unmarked boundary loops. 
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The relevant bare couplings in two dimensions are the bulk cosmological constant 
A2 and the boundary cosmological constants Zi and Zf. They are related to the 
renormalized continuum coupling constants A2, Aj, Aj by 



A2 = log2 + -A2a^, Zij = Aija. 



1 

It is straightforward to expand ([T7|) to lowest non-trivial order in the lattice spacing 
a and thus obtain the matrix elements {Af\H\Ai). We have 

(A;|f |A,) = analytic - log(A, + A^) - ^hz^^h±M + 0(0^), (19) 

where "analytic" refers to a constant term and to terms linear in Aj and Af. By 
an inverse Laplace transformation we can change variables from the boundary cos- 
mological constants Aj j to the conjugate (continuum) length variables Ljj of the 
boundaries. To first order in a, one finds 

-ioo ^''^ J -ioo ^''^ 

= l-6{U-Lf)-a(^^^+A^5{U~Lf), 

where we have ignored the analytic terms which lead to non-propagating terms of 
the form 6{Lf)6{Li). The appearance of such non-universal terms is familiar from 
the transfer matrix of 2d Euclidean gravity and can be accounted for by a careful 
treatment of the boundary conditions at L = ^ . 

In complete analogy with the definition given in footnote ^, the symmetry factor 
of a spatial S^-boundary of length L is given by Cl = 1/L and there is an orthogonal 
basis {\L)} of the continuum Hilbert space, obeying the normalization conditions 

(Li|L2) = ^(5(Li-L2), / dL\L)L{L\ = i. (21) 
J^i Jo 

We deduce that the Hamiltonian operator in the "L-representation" , acting on func- 
tions ip{L) = {L\ip), is 

H = -^L + A2L, (22) 

which is hermitian with respect to the measure LdL, as it should be. 

In principle we would like to use the same strategy to determine the quantum 
Hamiltonian H{A) and its spectrum also in the three-dimensional case. Of course 
we cannot be sure that this will lead to a simple differential operator as a function 
of the spatial volume, as was the case in c? = 2, (0). Even in two dimensions, there 
are models inspired by Lorentzian gravity where this is not the case [^6]. However, 
in the 2d cases studied so far it has always been possible to turn H into a local 
differential operator by transforming to a variable different from L. 
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Figure 2: A piece of a typical quadr angulation at t + 1/2. The three types of 
squares made from solid and dashed lines arise as sections of the (4,1)-, (1,4)- and 
(2,2)-building blocks. 

5 The matrix model correspondence 

We will now relate the Laplace-transformed one-step propagator (|15D to a matrix 
model. This will be done by showing that the two-dimensional configurations asso- 
ciated with a slicing at half-integer t of a sandwich geometry appear as terms in the 
perturbative expansion of a hermitian two-matrix model with AB AB-mteiaction. 
(A similar observation for the original Lorentzian 3d simplicial model was already 
made in [0.) 

Imagine a one-step geometry Q obtained by gluing the three types of building 
blocks of Fig. |l[ The intersection of this three-geometry with the spherical constant- 
time hypersurface at t -|- | can be visualized as a pattern of squares, whose edges 
correspond to the intersections of the time-like triangular faces of the 3d building 
blocks with this surface. Let us distinguish the two cases where the time-like triangle 
has its base either in the t- or the t + 1-quadrangulation by drawing its intersection 
at as either a solid or a dashed hne. Thus a (4,l)-pyramid gives rise to a square 
of solid edges, a (l,4)-pyramid to one of dashed edges and a (2,2)-tetrahedron to 
a square with alternating solid and dashed edges. The way in which these two- 
dimensional building blocks appear in the S^-quadrangulation at t+| is that they 
can be glued to each other only along pairs of edges of the same type. 

The quadrangulation at t + | can thus be viewed as a double-line graph of the 
kind illustrated in Fig. This type of graph is generated in the large-M limit by 
the two-matrix model 

Z{a,,a2,f3) = I dAdSe-^^*^(5^'+i5^-^^^-^s^-f^^^«), (23) 

where A and B are Hermitian M x M matrices. By expanding the non-Gaussian 
part of the exponential in powers of ai, 02 and f3 and performing the Gaussian 
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Figure 3: Matrix-model representation of the building blocks at t + 1/2. The gluing 
rules for the squares are determined by the Gaussian integrations, (AijAki) = SuSjk, 
{BijBki) = Sii6jk, and {AijBki) = 0. 



integral we are led by Wick's theorem to a successive gluing of the three kinds of 
squares described above if we make the index assignments as shown in Fig. |^. 

As usual, the logarithm of the partition function of the model, M^F(ai, a2, (3) = 
logZ(ai, a2, P), generates only connected quadrangulations, and taking the large-M 
limit will select those with S'^-topology. In principle F may be expanded in a power 
series in M~^, with higher-order contributions corresponding to quadrangulations 
of higher genera. Although we are mainly interested in the spherical limit, anything 
we say could be repeated for higher-genus surfaces, thus relating to 3d gravity on 
space-times with topology S*^^^ x [0, 1]. We can now write out the generating function 
F as an explicit power series, 

F(«i,a2,/?)= Yl Ar(iV4l,iVl4,iV22)«f«af^/3^^^ (24) 

N4,l,Nl4,N22 

where A/'(A^4i, A''i4, A''22) denotes the number of (connected) spherical quadrangula- 
tions described above, including symmetry factors. 

Comparing the form of ( P^ ) with the previous expressions (|r^) and 
making the identifications^ 

a, = e^-^-"% a2 = e^-^-"^ /3 = e-(5^+^), (25) 

one could be tempted to conclude that 

F(a^,a2,P) = {zf\f\zi). (26) 

However, this is not correct, since the number of configurations generated by the 
matrix model is strictly larger than those obtained from the Lorentzian gravity 
model, that is, 

■Ar(iV4i,iVi 4, iV22) > ^^{N,,, iVi4, iV22). (27) 

^ Adopting (p5|), the limit of vanishing boundary cosmological constants Zij corresponds to 
setting ai = a2- Unfortunately we cannot put ai =0.2 before having extracted the matrix elements 
(iV/|f|7V,) from F(ai,a2,/3). 
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Figure 4: Vertices of 0^-graphs dual to spatial quadrangulations. 



This difference is intimately connected to the fact that in the gravity case, the 
allowed configurations are not so much two-dimensional structures per se as sections 
of larger, three-dimensional objects whose three-dimensional manifold structure is 
encoded in the colouring (the dashed and solid lines) of the two-dimensional graph. 

The generalization inherent in the matrix model is best described in terms of the 
graphs dual (in a two-dimensional sense) to the quadrangulations. In this more con- 
ventional picture, the terms tr A'^, tr B'^ and tr ABAB are represented by four-valent 
vertices (placed at the centres of the squares of the original 2d quadrangulation) with 
four outgoing solid lines, four dashed lines or alternating solid-dashed lines (see Fig. 
I). The labels A and B are now associated with the dual solid and dashed edges 
connecting pairs of such vertices. 

In order to discuss the regularity conditions that must be satisfied by such a 
dual graph to qualify as (the dual of) a section of a Lorentzian three-geometry, we 
define a A-loop (a -B-loop) as a closed sequence of solid (dashed) dual edges with no 
further solid (dashed) dual links in its interior. This interior region (which has the 
form of a two-dimensional disc whose boundary is the loop) we call a A-domain (a 
5-domain). A dual graph coming from a 3d Lorentzian geometry then satisfies the 
following constraints (c.f. App. 2 of [p!8|): 

(1) The two subgraphs formed from only A-edges and only 5-edges must each be 
connected. 

(2) The two separate A- and -B-subgraphs can have neither tadpoles nor self- 
energy subdiagrams. This ensures that they are associated with regular 2d 
simplicial manifolds at the times t and t + 1 respectively. 

(3) The intersection of any pair of A- and 5-domains cannot be multiply con- 
nected. This implies that any pair of vertices, one at time t and one at time 
t -|- 1, of the original quadrangulation cannot be connected by more than one 
time-like link. 

(4) The (one-dimensional) intersection of a A-loop with a S-domain (and vice 
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versa) must be either empty or simply connected. This imphes that any three 
vertices of the original quadr angulation cannot belong to more than one time- 
like triangle. 

Thus we see that the dual graphs coming from the matrix model are considerably 
more general than those associated with 3d simplicial space-times. The generaliza- 
tions occur in several ways. The matrix-model graphs can have arbitrary numbers 
of disconnected spherical A- and i?-subgraphs (the only requirement being that the 
combined graph is spherical and connected) and each of these may contain (gen- 
eralized) self-energy diagrams and tadpoles. Furthermore, the A- and -B-loops can 
freely meander around each other, with arbitrary numbers of mutual intersections. 

We conjecture that the conditions (2)- (4) are not important in the sense that their 
implementation or otherwise will not affect the continuum properties of the model. 
We think that they constitute merely a 3d generalization of the universal behaviour 
already observed for 2d matrix models. In that case the inclusion of tadpole and 
self-energy subgraphs leaves the continuum limit unchanged (although it implies an 
enlargement of the configuration space from genuine 2d simplicial manifolds to more 
general 2d combinatorial complexes). 

We have checked numerically that the same happens in our 3d model. More 
precisely, we have performed simulations for the original simplicial model with tetra- 
hedral building blocks of |]T5|, TB| where one can formulate conditions completely 



analogous to (l)-(4) above. Dropping then the constraints (2)- (4), but keeping (1) 
we found that the key results of JlSl remained unaffected.0 

The status of condition (1) is different. Recall that in the original representation 
one obtains the geometry at time t (at time t + 1) from the quadrangulation at 
t + 1/2 by shrinking all dashed (solid) links to zero. The A- and 5-graphs of the 
dual picture at t -|- 1/2 are of course individually precisely the duals of these two- 
geometries at integer times. What does it imply for the spatial quadrangulation at 
t if the dual A-graph is disconnected? Using the same "shrinking-prescription" , it 
cannot separate into several pieces, but it does degenerate in the sense of forming a 
number of connected spherical graphs which touch each other pairwise only in single 
points. The resulting two-dimensional space is therefore no longer a manifold, but 
a branched tree of such spherical components. 

This phenomenon is illustrated in Fig. ^. The pictures at the top show two quad- 
rangulations of a two-sphere at time t-|-l/2 made of the three types of square-shaped 
building blocks. As two-dimensional configurations they look rather innocent and 
regular. However, they are pathological from a three-dimensional point of view. 

^^However, one interesting change did occur: in ||l8| we observed a first-order phase transition for 
large k (there called feg) to a phase where successive spatial slices decouple. We viewed the presence 
of this phase as a discretization artifact irrelevant for the continuum physics. This interpretation 
is corroborated by our recent simulations without the constraints (2)- (4), where the large- A: phase 
is simply absent. It implies that the generalized model gives in some sense a better representation 
of continuum physics. This is reminiscent of 2d simplicial quantum gravity, where the continuum 
limit is approached much faster if self-energy and tadpole graphs are allowed. 



15 



(a) 



(b) 



Figure 5: Examples of matrix-model configurations at t + 1/2 which are not allowed 
in the original Lorentzian gravity model and which result in geometries with worm- 
holes at time t. Shrinking the dashed links to zero, one obtains the two-geometries 
at the bottom. The thick dashed lines in the quadr angulations at the top are con- 
tracted to touching points or to points along one-dimensional wormholes. 

Drawing the dual graphs, one finds in both cases that the A-subgraph, made of 
solid dual edges, consists of two components. Each component is dual to a quad- 
rangulation of a two-sphere by four (solid-line) squares. When shrinking away the 
dashed links of the original quadrangulation, the "necks" (indicated by the thick 
dashed circular lines) between the two spheres are gradually pinched to points. If 
there is just one neck, the two spheres will touch in a point (Fig. |^a). If there 
are several concentric necks, with annuli of dashed-solid squares (corresponding to 
closed rings of (2,2)-tetrahedra) in between, this process generates one- dimensional 
"wormholes" , as illustrated in Fig. ||b. 

Similar so-called "touching" -interactions have been studied in the context of two- 
dimensional matrix models |2^, ^ , and are closely related to these wormholes (see 
Sec. ^ for a detailed discussion). Although configurations of this type are already 
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present in ordinary matrix models of Euclidean 2d gravity, the explicit introduction 
of such an interaction term in the action (and an associated coupling constant) al- 
lows us to increase their weight and thereby control their abundance. Depending 
on the value of the coupling constant, one can obtain either the ordinary universal 
behaviour of two-dimensional Euclidean gravity, or a modified critical behaviour, 
described in more detail below. We will see that in our 3d-gravitational reinterpre- 
tation of the matrix model the presence of wormholes and touching-interactions is 
governed by the value of the bare gravitational coupling constant k. 

Let us finally note that the above discussion about degenerate geometries could 
be repeated verbatim for the cubic matrix model with partition function 

Z{a^,a2,f3) = J dAd5e-^*^(5^'+i^^-^^^-¥^^-f^^^«), (28) 

which is associated with the original 3d simplicial gravity model described in [|1^, |1^] 
(with (3,1)- and (l,3)-tetrahedra instead of (4,1)- and (l,4)-pyramids). The only 
reason why we prefer to use the quartic model defined by (|2^) is that it has been 
solved for ai = a2 We would expect from universality arguments that the 

models given by (p3| ) and ( pHD lead to the same continuum physics. 



6 Reinterpreting the matrix model 

We have seen in the previous section that 3d Lorentzian "sandwich" geometries 
can be put into correspondence with a subclass of 2d graphs generated by a quartic 
matrix model with Ai^Ai^-interaction. We have also argued that this generalization 
is potentially relevant, in the sense of the two models having a different phase 
structure. 

The phase structure of the three-dimensional Lorentzian model has been inves- 
tigated numerically in For given gravitational coupling k, there is a critical 
value Xc{k) which is to be approached from above, that is, from the region A > Ac 
of the cosmological constant where the partition function converges. In the orig- 
inal simulations we also observed a first-order transition for large k when moving 
along the critical line Xc{k). As already mentioned in footnote |l^, this transition 
is a discretization artifact which disappears when one relaxes some of the mani- 
fold constraints. What emerges as the phase structure for 3d Lorentzian gravity is 
therefore simply that of a single phase (at least in the range of coupling constants 
accessible to our computer simulations) where the taking of the continuum limit co- 
incides with the tuning of the cosmological constant to its critical value, much as in 
two-dimensional models of quantum gravity. The role of the gravitational coupling 
k is merely to set an overall scale for the system, without affecting its continuum 
properties, as we have argued in 



What we want to explore presently is the physical interpretation in terms of 
three-geometry of the generalized model (where multiply-connected A- and 5-graphs 
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Figure 6: The phase diagram of the ABAB-model, according to with the critical 
hne and the critical point at a = /? = (corresponding to k = 1). One way of 
approaching the critical line is through fine-tuning of the cosmological constant A 
along lines of constant k, shown as dashed curves. To end up at the critical point, 
one should move along the line of constant kc = — | In j3c — 0.85. 



are allowed to occur as the spatial sections at half-integer times) and of its phase 
structure. This is of course the two- matrix model defined by eq. (pSj). Let us consider 
first the case of vanishing boundary cosmological constants, Zij = 0. Although we 
will not be able to construct the Hamiltonian or the propagator (c.f. footnote ^), 
we can nevertheless discuss the phase diagram of the model. We will comment later 
on the general matrix model with Zi, Zf ^ 0. 

Physically this choice of couplings implies that we are studying the geometric 
fluctuations between two successive spatial slices with free boundary conditions. 



As was clear from the computer simulations in |18[, the behaviour of this one-step 
system determines the phase structure of the discretized theory. Other than that, 
the choice Zij = has the great advantage that the corresponding matrix model has 
already been solved. Details of the solution can be found in the paper of V. Kazakov 



and P. Zinn- Justin ^0|. Their analysis of the matrix model with ai = a2 = a can 
be summarized as follows. Let us fix a ratio k = a//3 and gradually increase the 
values of a and (3 away from zero. The large-M limit of the matrix model is defined 
for sufficiently small a and (3, and the model becomes critical at a point (/5c(/t), 
ac{n) =k/3c(k)), giving rise to a critical line in the a)-plane (see Fig. Along 
it one finds two phases, separated by a second-order phase transition at k = 1, 
with ac{l) =/3c(l)- In the context of two-dimensional gravity, they were given the 
following interpretation [^: the phase of k > 1 with ad^) > /?c(^) can be viewed as 
ordinary two-dimensional Euclidean gravity, with central charge c = 0. In this phase 
the relation between the length I of a "typical" loop of A-links (or 5-links) and the 
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area A of the two-dimensional surface enclosed by the loop is ~ A. By contrast, 
in the other phase (which also has c = 0), the scaling is anomalous, /^/^ ~ A. 

One can understand these two situations by simply looking at the two limits 
(3 = and a = 0. By setting (3 = we switch off the Ai^Ai^-interaction, leading 
to a decoupled system of two (p^ one-matrix models, whose critical behaviour is 
individually that of an ordinary 2d Euclidean gravity system. Setting a = 0, we 
obtain a matrix model whose only interaction comes from the term tiABAB. This 
model was first solved in by mapping it to the so-called dense loop phase of the 
0(l)-matrix model on random graphs.[^ 



6.1 Intermezzo: touching- interactions 

We will now reinterpret the results described above in the context of three-dimensio- 
nal Lorentzian gravity. Recall our earlier discussion in Sec. ^ where we showed that 
from a three-dimensional point of view a generic feature of the matrix-model config- 
urations is the presence of touching points and wormholes in the associated spatial 
slices. The matrix model therefore describes transitions between two-geometries 
with a tree-hke structure, consisting of regular quadrangulations of two-spheres, 
pairwise connected by one-dimensional wormholes of length / > (Fig. ^). 

Let us for a moment revert back to a pure matrix-model language. For the sake 
of definiteness, we will concentrate on a quadrangulation at time t made from solid 
squares. We will now show that the effective weight associated with the mutual 
touching points of the two-spheres making up the tree-configuration is obtained by 
integrating out the S-matrices. While this cannot be done explicitly as long as the 
tT term is present, it is instructive to do the integration for the general matrix 
model (|23|) with a2 = 0.|3 In this case, the action is Gaussian in the i?-matrix and 
we obtain 

f Af ( -itr(yl2+B2)+^trA4 + i/3 trAB^B ) 

Z(ai,a2 = 0,/?) = dAdBe \ ' ' ' J 

= dA (^det{I - PA^ A)j e V ^ ' J 

f A/( -itrA2 + ^^trA4 + -i^ V?° , ll^tr AHr aA 

= dA e \ ' * ' /, (29) 

which now describes a particular one-matrix model. Terms like tr A'^tr A'^ are usu- 
ally referred to as "touching-interactions". How can their effect be visualized in 

The 0(l)-model on random graphs is known to correspond to c = quantum gravity, but 
with two distinct phases |3l|, not unhke those of the ABAB-mode\. In a continuum 
interpretation, these result from assigning different boundary operators to the theory. The dilute 
loop phase corresponds to ordinary 2d Euclidean quantum gravity where boundary lengths scale 
according to P ^ A, while in the dense phase one has l^/^ r^A^,^. 

^^For 02 7^ there will be other touching interactions beyond the ones described below. 
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Figure 7: A typical configuration generated by the matrix model ( p3D for 02 = 0. By 
integrating out the i?-matrices one obtains touching terms tr A'^tr A'' (here, k = A 
and 6). Their effect is to identify the boundary links of the solid-square components 
pairwise across the "gaps" formed by the solid-dashed squares, as described in the 
text. 



geometric terms? Firstly, in the same way as ti A^ represents a square of solid 
edges, a term tr A^ in the action can be thought of as a fc-gon. Expanding the 
action and performing Wick-contractions corresponds to gluing such polygons to- 
gether by identifying their links pairwise. Imagine that during this process we had 
constructed a surface glued from various polygons, with a boundary consisting of k 
links. By performing the appropriate Wick-contractions, we could now glue a fc-gon 
(represented by tr A'') to this boundary and close off the surface. 

By extension, we can determine the effect of a term like trA^trA^. Imagine 
a pair of surfaces made from polygons, both with a boundary of length k. By 
contracting with tr A'^tr A'' , we can close off the two surfaces. Because this happens 
simultaneously for both of the surfaces, we associate by this process a fc-gon of the 
first surface with another /c-gon of the second surface, which we may think of as a 
touching point between the two geometries (for example, located at the centres of 
the fc-gons). 

This is illustrated by Fig. 0, which shows a typical quadrangulation at t + 1 /2 
appearing in the matrix model with a2 = (so no dashed squares are present), 
before performing the S-integration. The shaded areas represent components made 
of solid squares, where the outer boundary should be thought of as a single point, 
since the entire geometry is spherical. Integrating out the S-matrices has precisely 
the effect of generating terms tr A'^ti A'' and identifying the solid boundary links 
pairwise (each link with its partner lying opposite in a solid-dashed-solid-dashed 
square). Within the one-matrix model it is of course a matter of convention how 
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one wants to think of the pairwise associations introduced by the terms tiA^trA^, 
whether as microscopic touching points between the two closed-off surfaces or simply 
as regular gluings of one boundary of length k to another. However, it is natural 
to think of them as creating a connection between the surfaces. There is a factor 
1/M associated with each touching interaction, and it is well-known that to leading 
order in M one obtains through "touching and connecting" only two-geometries 
with an overall spherical topology, and no spaces of higher genus. Collapsing all the 
genuinely two-dimensional parts of this two-geometry to one- dimensional lines, we 
would obtain a branched polymer, i.e. a one- dimensional tree structure without any 
closed loops. 

What is our motivation for this excursion into touching-interactions? The general 
ABAB-model does not in itself suggest an immediate physical interpretation in 
terms of spatial or space-time geometry. However, with the benefit of hindsight, 
one could have taken the presence of the limiting case /3 = (where one has two 
independent copies of Euclidean 2d gravity) as an indication that switching on the f3- 
coupling will lead to a theory of two interacting two-geometries. The situation where 
such an interaction occurs naturally is of course when the two-geometries appear as 
neighbouring embedded spatial slices within a 3d space-time, with their interaction 
dictated by the 3d Einstein action. Given the three-dimensional interpretation, it is 
then natural to determine the effect of the neighbouring two-geometry by integrating 
out the corresponding degrees of freedom in the partition function (for instance, the 
i?-matrices) . As we have seen, this effectively introduces touching-interactions in a 
2d slice which are dictated by (3, where (3 in turn depends on the coupling constants 
A and k of three-dimensional quantum gravity. 

Another related point that follows from the three-dimensional interpretation is 
that it gives us a precise geometric way of obtaining a two-geometry at time t (or 
t -|- 1) from a general matrix-model configuration at time t -|- 1/2 by the shrinking- 
prescription already mentioned in Sec. |^. This may be visualized as a continuous 
process where the lengths of either the dashed or the solid links are gradually shrunk 
to zero. Although this prescription was originally invented to describe the well- 
behaved geometries of 3d Lorentzian quantum gravity, it works just as well for the 
more general configurations generated by the matrix model. We have already seen 
explicit examples in Sec. |^, where this led to the creation of spatial wormholes (Fig. 

I). 

What is important to notice is that from a three-dimensional point of view, 
such wormholes and the touching-interactions described in this subsection are really 
two sides of the same coin. That is, taking a general configuration of the ABAB- 
model which is not allowed in the original Lorentzian model, one will find both types 
of degeneracy when applying the shrinking-prescription. We illustrate this with a 
simple matrix- model graph in Fig. |^. We start at time t + 1/2 from an extended 
version of Fig. ||b, where the shaded areas indicate regular dashed quadrangulations. 
Since these will all vanish when we shrink away the dashed edges, we obtain the 
same wormhole geometry at time t as before. On the other hand, if we now shrink 
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t+1/2 

Figure 8: A configuration of the ABAB-matiix model, giving rise to a two-geometry 
with a spatial wormhole at time t and a two-geometry with "touching-interaction" 
at time t+1. Viewed as part of a Lorentzian space-time, this configuration is 
degenerate. 

the solid edges to zero, we are left with two pieces of dashed quadrangulations, each 
with a boundary of dashed links of length four (the thick dashed lines in Fig. |]) . We 
then must identify the boundary links pairwise, as indicated by the arrows. This 
is clearly the same mechanism that was at play in the touching-interactions, now 
generalized to a situation with 02 7^ 0. Although the resulting configuration at time 
t + 1 is (unlike the wormhole at t) perfectly well-behaved as a two-geometry, it is 
not acceptable from our original three-dimensional point of view because of the way 
it is "embedded" in the three-dimensional space-time. We will in the following not 
distinguish between those two types of geometric degeneracy and simply refer to 
them as "wormholes" 

6.2 Phase structure in the presence of touching-interactions 

Having established the relation with one-matrix models containing (generalized) 
touching-interactions, let us try to understand what we can learn from the known 
phase structure of these models. There are three phases when touching-interactions 
are added to the usual tr A^-interaction. In phase I the touching-interaction plays 
no role, and the critical behaviour is determined by the tr A^-interaction. A typical 
configuration (with probability one in the scaling limit) is a connected, spherical 
quadrangulation of the type that dominates Euclidean 2d gravity. 

^•^This fits with the terminology used in matrix models, where the touching points brought about 
by the touching-interactions are also sometimes thought of as "wormholes" . 
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In phase II, the touching- and the tr ^^-interactions coexist and compete. A 
typical two-geometry consists of many spherical 2d "baby" -universes of all sizes, 
which are connected to each other at touching points, thus forming a "blown-up" 
branched polymeij^. As was first realized in |3^, the critical behaviour results 



from an interplay between the properties Euclidean gravity and those of branched 
polymers (see also |p5|). 



Quite surprisingly, all critical aspects of this phase can be understood from 
continuum quantum Liouville theory (which is usually thought to describe only 
phase I) p6| , [37| , p8[ . Let us review this briefly. Consider a conformal field theory 
in dimension two with a fixed background metric g, and let $ be a spinless primary 
field with scaling dimension A^^\ Then the "one-point" function for $ scales as 

Fi'\A) = ( /d^e vW)m) - (30) 

\>'' / CFT 



where A = Jj. d^^ \/ g{0 is the area of the underlying manifold E, and the ex- 
pectation value is taken in the conformal field theory. Coupling this theory to 2d 
Euclidean quantum gravity, the metric g will be allowed to fluctuate. We can decom- 
pose a general metric g into g = e'^g, where 0(^) is the conformal field. In conformal 
gauge, the integration over the fluctuating metric becomes an integration over the 
conformal factor (j), weighted by the Liouville action >S'l[<^], while the spatial integral 
appearing in the one-point function is changed according to 

I d'^^/mm I d2e\M)e'^**(«)$(0. (31) 

The dressing exponent /5$ is determined by requiring that the "dressed" operator 
g/3*0,|) j^ave conformal dimension dimension (1,1), so that it can be integrated over 
the two-dimensional surface without breaking conformal invariance. This leads to 
the relation 

(0) . [25 



/3$(/3* + g) = -2 + 2A^^^ Q = y__. (32) 

The special case of the unit operator $ = 1 corresponds to A^^"* =0, with dressing 
exponent (3i. The one-point function for the gravity-coupled theory is defined by 



QG 

-.^WrdV^/^e^*'^cI)\„.__ 

, (33) 



/ 1)0 eS^^-l^] 6(1 d'^^e^^^-A) (/ d^^ v^e^* ^ <l>>^^^ 



ZiA) 

where Z{A) is the partition function with fixed area A (enforced by including the 
same (5-function that appears in the numerator of (|33[)). By requiring that F^{A) be 



^Further information about the properties of branched polymers can be found in |32 
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independent of the fiducial background metric g, a change of integration variables 
— >• + ^ log A leads to 

F^{A) = A"-^'^' (34) 

In analogy with (|30|) , we now define the critical exponent A$ of $ for the gravity- 
coupled theory by 

F^{A)^A^-^^. (35) 



Combining this with relation (p4f), one derives 



A. = 1 - ^. (36) 
Pi 



For our present purposes it is important to note that eq. (|3^) has two solutions, 
namely, 

_ -(25-c)±^(25-c)(l-c + 24Af^ 

In conventional quantum Liouville theory one works with (3^^^ because this choice 
ensures that in the "classical" limit c = — oo (classical in the sense that the fluctua- 
tions of the Liouville field in (^) are completely surpressed) the scaling dimensions 
of the model without gravity coupling are recovered, 

A$ ^ for c ^ -cx), (38) 



as can easily be verified from formulas (|37| ) and (0). 

It is a remarkable fact that all critical exponents which can be calculated in 
phase II of the matrix model with touching-interactions are obtained by replacing 
/^i^^ ~^ changing the scaling relations according to 

F^^\A) ~ A^-'IP'^"' — F^^\A) ~ A^-'l^'^"\ (39) 

The implications of this prescription can be illustrated by considering the simplest 
boundary operator, 

t = I ds^I{s), (40) 

JdT. 

which measures the length of the boundary of the two-dimensional manifold E as a 
function of the induced metric e(s) on 9S. Unlike in matrix models, it is somewhat 
awkward to introduce boundary operators like (^Oj) in Liouville theory. Nevertheless, 
C. can be treated along similar lines as the operator $ in (|3lD above. The expression 
one obtains for the gravitationally dressed version of the length operator is 

dsv/i(^e^^"'^(^), /?W = \(3^^\ (41) 
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From the point of view of conformal field theory it is rather surprising that the 
"naive" value (3^ = \j3i (i.e. the boundary scales like the square root of the bulk) 
is correct. However, as can be read off from (0), this does lead to the canonical 
scaling 



(£(+)(^))Q^^^/3rv/^i=^i/2 (42) 

already known from matrix-model considerations. 

To determine the scaling relevant in phase II, we should now replace the operator 
It turns out that the correct substitution for this particular boundary operator 
is given by 



This leads to the scaling behaviour 



/^r — /3r = 2/5r- (43) 



{1^''\A))qg-A^'^''/^'^' =A''\ (44) 

implying that the boundary scales anomalously in terms of the area. This anomalous 
scaling is clearly a reflection of the additional fractal structure introduced by the 
touching- interactions. Note that it coincides with the spectral dimension of branched 
polymers 



The scaling (^) coincides with the scaling behaviour observed in the ABAB- 
matrix model in the phase with k < 1, and mentioned already earlier in Sec. ^ This 
corroborates our interpretation of this model in terms of interacting two-geometries. 
Furthermore, it agrees with the relation between loop length and area found in the 



dense phase of the 0(l)-model, where the replacement (|^) was first noticed [p^ . 

Lastly, the critical behaviour in phase III of the matrix model with touching- 
interactions is characterized by a complete dominance of branched polymers. The 
size of the individual spherical components making up the 2d universe never exceeds 
the cutoff scale. This phase does not seem to have an analogue in the AB AB-m.a.tiix 
model. 



6.3 Summary of Sec. |^ 

Let us summarize what we have learned about the three-dimensional geometric 
interpretation of the ABAB-model. To this end, we will first translate its phase 
diagram (for ai = to the A;-A-plane of the gravitational couplings. We deduce 
from (^) that 

— = e3'= and — ~A; = -log— -, 45 

Gn 3 (3c{k) 

where the second relation holds along the critical line (the points where a non- 
trivial continuum limit exists). The shaded areas in the two phase diagrams of Fig. 
^ indicate the region in which the partition function converges. The dashed curves 
in both diagrams show an approach to the critical line along a curve of constant k. 
Moving from left to right along the critical curve of the left phase diagram (so that 
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Figure 9: Translating the phase diagram of the AB AB-matrix model in terms of 
the original coupling constants a and (3 (left) to the couplings k and A associated 
with its interpretation in terms of 3d Lorentzian gravity (right). 

f3 increases and n decreases) translates into a motion from right to left along the 
critical line of the "gravitational" phase diagram (that is, towards a smaller k or a 
larger bare Newton's constant Gat). 

In terms of space-time geometry, at the point Gn = (corresponding to (3 = 0) 
subsequent spatial slices are completely decoupled. No information can be propa- 
gated in time and the system has no interesting three-dimensional properties. In- 
creasing Gn away from zero, we enter the phase k > 1 of the matrix model. Interac- 
tions between neighbouring spatial slices of the 3d universe become possible|3 and 
the geometrical properties of a typical spatial slice at integer-t resemble those gen- 
erated in 2d Euclidean quantum gravity. It is very suggestive to identify this phase 



with that of Lorentzian gravity defined in |T^, 0, since the simulations of [|1^, |1^] 
revealed both the existence of extended three-dimensional space-times (indicative 
of correlations in time-direction), and evidence that the Hausdorff dimension of the 
spatial slices is dn = 4, in agreement with that of 2d Euclidean gravity. 

Apparently, the phase with k < 1 cannot be realized in the Lorentzian manifold 
model because it does not allow for the creation of "wormholes"[^. By contrast, 
no such restriction exists in the ABAB-ma.tTix model, and by gradually increasing 
Newton's constant Gn- such wormhole configurations become energetically more 
favoured. Looking at spatial slices in the phase k < 1, a typical two-geometry will 
consist of many smaller universes connected to each other by cutoff-size wormholes. 



^■^We cannot be more specific about the nature of the interactions before having analyzed the 
transfer matrix of the model in more detail. 

^®At least, wc have so far not seen any evidence of this phase or the associated second-order 



transition in the computer simulations. 
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7 The general ABAB-matrix model 



While the phase structure of the matrix model of (generalized) 3d Lorentzian gravity 
can be understood from the KZ-solution for the special case ai = of the ABAB- 
matrix model, the construction of the quantum Hamiltonian requires that we perturb 
away from ai = 0:2. It is not difficult to generalize the ansatz of [EH] to this situation, 



using the techniques developed in ||41[]. One can write down a set of singular integral 
equations which generalize those solved in |]2D|. Furthermore, using results from 
the general theory of singular integral equations one can prove that the solution is 
unique. 

Although we have at this stage nothing much to say about the explicit solution 
for general ai and 02, it may be worth pointing out a connection with yet another 



2d statistical model. As mentioned in |2^ and studied in more detail in |i2], the 



AB AB-raatiix model with ai =ol2 can be mapped to an eight- vertex model defined 
on a random four-valent lattice (in the dual picture). The map is simply given by 
X = A + iB. Using the same transformation, one can map the general ABAB-model 
to a sixteen-vertex model with partition function 



Z{ai,a2,P)=jdXdX^e 



Ntr 



where the constants b, c, d, and / are given by 

ai + a2 + 2/5 ai + 02 - 2/? ai - 02 

b = , c = a = , f = . 46 

16 16 16 ^ ^ 

Since the eight-vertex model on a regular square lattice has already been solved, it 
may not be very surprising that the restricted ABAB-model can be solved too. This 
would be in line with the general observation that it is often simpler to solve matter 
models on random dynamical lattices (associated with 2d quantum gravity) rather 
than on regular ones. By the same token, one should maybe not be too discouraged 
by the fact that the general sixteen-vertex model has not been solved on a regular 
lattice. This is also not what is needed here, because of the particular form of the 
parameters (|^). 



8 Taking the continuum limit 

We have already described in some detail in Subsection ^]3| how the phase structure 
of the matrix model with ABAB -inter action can be reinterpreted from the point of 
view of three-dimensional geometry, and how this fits into our previous investigations 
of three-dimensional Lorentzian quantum gravity. One further point we would like 
to make concerns the nature of the continuum limit in the gravitational model, 
and what the newly established relation with the matrix model can add to our 
understanding of it. 



27 



In dynamically triangulated models of quantum gravity one usually performs 
the continuum limit by fixing the bare (inverse) gravitational coupling constant k 
and fine-tuning the bare cosmological constant A to its critical value Xc{k). Since 
the space-time volume is conjugate to the cosmological constant in the action, such 
a fine-tuning corresponds to taking the lattice volume Ntot = Ni4^ + A^4i + ^N22 to 
infinity, and can be viewed as an additive renormalization 

A = Ae(^) + AA, AA = Asa", (47) 

of the cosmological constant, where A A should be related to the continuum cosmo- 
logical coupling constant A3 by a suitable scaling. 

The analogue of this procedure in a lattice field theory, for example, the Ising 
model defined on a hypercubic A^'^-lattice, would be to take the infinite-lattice limit. 
However, this may not lead to a continuum limit. In the case of the Ising model, 
in order to obtain such a limit, the coupling constant (3 (the inverse temperature) 
must always be fine-tuned to a critical value, even if the lattice volume was infinite. 
By analogy, one might therefore expect that also in simplicial gravity models, the 
gravitational coupling k had to be fine-tuned to a critical value kc to arrive at an 
interesting continuum theory. 

However, this is not the only possible scenario. Consider another hypercubic 
lattice model in d dimensions, one with a lattice scalar field 0„ and the simple 
Gaussian action 

d 

S[<P] = Y.Y.^<Pn+e.-<Pn?. (48) 

n 1=1 

For this model, one automatically obtains the continuum Gaussian field theory by 
taking the limit of infinite lattice volume and rescaling the lattice spacing to zero. 
The same situation is encountered in two-dimensional simplical quantum gravity, 
both for Lorentzian and Euclidean signature, when formulated as a sum over tri- 
angulations with geodesic edge lengths determined by a single lattice spacing a. 
Fine-tuning the cosmological coupling constant A^^-* to its critical value according 
to 

=A(2)+AA(2), AA(2)=A2a2, (49) 

or, equivalently, taking the infinite-volume limit and scaling a to zero, leads auto- 
matically to a continuum theory of 2d quantum gravity. In (^), A2 denotes the 
continuum two-dimensional cosmological constant, which has already appeared in 
our earlier discussion of Lorentzian 2d gravity (see eq. (|TB|)). 

It is easy to see that the taking of the continuum limit in the two-matrix model 
with AB AB-interaction falls into this latter category. According to EH], the con- 



tinuum limit is obtained by fixing the ratio K = a/ (3 (corresponding to straight lines 
through the origin in Fig. ^), and increasing a and (3 until the critical point is 
reached. This is a natural procedure to adopt, since the ratio a/ [3 appears in the 
equations which determine the solution of the model. However, one could equally 
well approach the critical curve ac(fi;) = n(3c{K) along any family of curves in the 
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{(3, a)-plane which intersect this hne. In particular, one could choose curves of con- 
stant k, which according to ( ^5] ) are defined by a = k(3'^ (c.f. Fig. The curve for 
a given value k intersects the critical line in the point {f]c{i^{k)), ac{K{k))). Moving 
along the constant-A; curve to the critical point is of course tantamount to fine-tuning 
the cosmological constant, described by eq. (^Tf). 

This observation lends further support to our interpretation of the role of the 
gravitational coupling constant k in 3d Lorentzian quantum gravity. Based on the 
results of our computer simulations we suggested in that the continuum limit of 
this theory does not involve a renormalization of k. Rather, k should act as a finite 
overall scale parameter. Now, as we have seen in Sec. it is natural to identify 
the Lorentzian model of fl^, |l^ with the k > 1 phase of the ABAB-matiix model. 
Forgetting for a moment about the new three-dimensional interpretation we have 
given it, the matrix model by itself is simply a two-dimensional statistical system 
with two phases and a transition point in between. In one of the phases, the scaling 
behaviour of the subgraphs of A- or i?-type corresponds to that found in ordinary 
2d Euclidean quantum gravity with c = 0. The other phase corresponds also to 
c = quantum gravity, but with the "unconventional" branch of the KPZ equation 
chosen in the way described in subsection p.2\ Changing k (that is, changing k), 
but staying within the same phase, the continuum limit is of course identical, which 
means that the theories can differ at most by the choice of an overall length scale. 



9 Outlook 

As we have already mentioned earlier, extracting further physical information about 
3d Lorentzian gravity, with or without wormholes, requires a more detailed under- 
standing of the transfer matrix (or, equivalently, of the quantum Hamiltonian) . This 
should help us to determine the appropriate renormalization of the boundary cos- 
mological constants in the continuum limit. We have seen in Sec. ^ that in the 
two-dimensional Lorentzian model, where the transfer matrix, the Hamiltonian and 
their spectra are known explicitly, the analogous scaling is canonical. 

Another interesting question is how the presence of the extended three-dimen- 
sional structures observed in simulations of the Lorentzian model is reflected in the 
mathematical properties of the transfer matrix. We found in |TB[ that when the 
extension in time-direction of the simulated three-geometries exceeds A3 - with 
the renormalized cosmological constant A3 setting an effective volume scale of the 
system - a universe with a flnite spatial volume exists only for a time duration 
T ~ A^^/^, whereas during the remaining time the size of the (compact) spatial 
directions remains at the cutoff scale. 

Our interpretation of the ABAB-mathx model raises the interesting issue of 
whether the phase transition of this model can be given a meaning in terms of 
three-dimensional geometry too. The intriguing fact that a two-dimensional central 
charge c = 1 is associated with this point brings to mind attempts of associating 
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2d conformal field theories with two-dimensional (boundary) surfaces in theories of 
3d gravity. Examples are given by the recent entropy calculations in anti-de Sitter 
43| and de Sitter space WM. We cannot directly compare with any of these results, 



because our set-up is quite different. We are working in a geometric and not a 
gauge-theoretic Chern-Simons formulation, our boundary surfaces are compact and 
space-like, and our (bare) cosmological coupling constant is by necessity positive. 
Nevertheless, the ABAB-matrix model is suggestive of the presence of additional 
symmetries at the phase transition point. 

As for the topology of our model, we have restricted ourselves to discussing space- 
times ^2 X [0,1], since the solution of the ABAB-mathx model as presented in 



is valid only for spherical topology. In order to make a comparison with other ap- 
proaches to (2-1-1) quantum gravity, it would be desirable to consider also surfaces of 
higher genus g, where the physical configuration space is described by a finite num- 
ber of Teichmiiller parameters. Results obtained so far for the matrix model suggest 
that its phase structure does not change for higher genus; the values of the critical 
coupling constants are independent of g, and most critical exponents are unaltered 
since they refer only to the short- distance behaviour of the model. Nevertheless it 
will be difficult to compare with a canonical reduced phase space quantization, say, 
because in the matrix model all modular parameters appear integrated over for a 
given genus. It is possible one could find a way to probe the individual Teichmiiller 
parameters, although it seems already quite complicated to perform the necessary 
large-M expansion of the AB AB-matiix model using the character expansion of 



0. 



We have argued in this article that under certain assumptions the two-matrix 
model with ABAB-intei action describes three-dimensional Lorentzian quantum gra- 
vity. In fact, not only does it describe the "regular" Lorentzian quantum gravity 
formulated in |T^, but a more general theory where space can split into many 
components connected by "wormholes" , resulting in a tree structure of spatial uni- 
verses forming a branched polymer. Starting at vanishing bare Newton's constant 
Gat, by increasing Gjy gradually we first find a phase of "weak gravity" which seems 
to coincide with the phase of regular Lorentzian gravity seen previously in computer 
simulations. Eventually, we meet a second-order transition point beyond which lies 
a phase of "strong gravity" and large Gn- From our original point of view of using 
causality as an effective regulator of quantum geometry, it is of course unclear to 
what extent the space-times with an abundance of spatial wormholes found in the 
strong-gravity phase are still acceptable in the path integral. In order to decide this 
question, we need to get a better understanding of the genuinely t/iree-dimensional 
properties of this system of quantum geometry. 
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Appendix 

In this appendix we will calculate the Regge version of the Euclidean gravitational 
action (^, which we use in Sec. 0. It should be remembered that our geometries 
are Wick- rotated versions of discrete Lorentzian (2+l)-dimensional space-times |1^ . 



For a change, we will do the angle calculations for Euclidean signature, and rotate 
back afterwards. Let the ratio of the squared lengths of time- and space-like links 
be given by a constant^] /? > |, such that /j = pi^. The dihedral angles of (4,1)- 
and (l,4)-pyramids around space- and time-like links will be denoted by 6s and 6t 
respectively. Calling the corresponding dihedral angles of the (2,2)-tetrahedra (ps 
and (pt, one derives the relations 

cos^ 9s = - cos = 4^^_ ]^ ' (ps = T^-'^ds, (pt = 7^-0t. (50) 

The numbers of the three different types of building blocks filling out the space- 
time between t and t + 1 are A^i4(t), iV4i(t) and N22{t). From this we compute the 
number of time- like links between t and t+1 as N^^{t) = Ni4{t)+N4i{t)+N22{t)+2. 
The number of spatial links contained in the spatial slice at t is 2A^4i(t), and each 
spatial link at t belongs to two (4,l)-pyramids. Collecting all this information and 
using the well-known expressions for Regge curvature in terms of dihedral angles 
including boundary terms EO], the curvature contribution to the total action 



is found to be 

- / d^x^R{x) + [ d^x^/hK{x) — > (51) 

2 Jm JdM 

^ a(47r + [(27r - 46^)^+ (27r - 4^,)] {N,,{t) + N^,{t) - N22{t))) , 



where the term in square brackets is always positive. The total volume contributing 
to the cosmological term in the action is given by 



l/^Va — " ^3 ' E (^4l(t) + iVl4(t) + ^iV22(t)). (52) 



^^The constants a and /? used in this appendix have nothing to do with the couphngs a and /? 
appearing in the main text. 
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Putting in the appropriate (bare) coupling constants, this leads to a total discretized 
action 

S[k, A] = Cot - ^iV4i + Nu - N22) + A(iV4i + Nu + ^^2), (53) 

where now the total numbers of the different building blocks appear and where the 
constants are given by 



/ ^ ^ i y ~ \ 

A ^ - + 2a//3 arccos — + arcsin — t— 

AtiGn\ ^ 4/3-1 (3-\ 



In order to obtain the Lorentzian version of the action (|53|) (up to an imaginary 
factor — i), we need to replace /? —a such that = —otll, in accordance with 
the notation in |15|, |T6[. To get the correct Lorentzian action, one then inverts the 



prescription given in [T^, |T6| by continuing the square-root expressions according to 



a/5 and \J l3 — \ i\Ja + |. 
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